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0> ' We analyze the generation of seed magnetic fields during de Sitter inflation consid- 

o : 

■ ering a non-invariant conformal term in the electromagnetic Lagrangian of the form 
^ ' 

5h ! —\l{(j))Fp^pF^^^^ where is a pseudoscalar function of a non-trivial background 
P-i. 

I field (f). In particular, we consider a toy model, that could be realized owing to the 



coupling between the photon and either a (tachyonic) massive pseudoscalar field and 



. a massless pseudoscalar field non-minimally coupled to gravity, where / follows a 

Q \ simple power-law behavior I{k,r]) = g/{—kri)^ during infiation, while it is negligibly 
H 

• small subsequently. Here, g is a positive dimensionless constant, k the wavenumber, 



rj the conformal time, and /3 a real positive number. We find that only when (3 = 1 
and 0.1 ^5^2 astrophysically interesting fields can be produced as excitation of 



> 

jC^ I the vacuum, and that they are maximally helical 
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I. INTRODUCTION 



It is not excluded that large-sca 
galaxies have a primordial origin [l, 



€, microgauss magnetic fields detected in any type of 
2|. If so, they have probably been generated during an 



infiationary epoch of the universe, since in this case their correlation length can be as large 



as the galactic one. (For possible generating mechanisms see, e.g., Ref. 0, Q, B, B-) An 
overwhelming proof of "primordial origin" of galactic fields could come from the observation 
of peculiar imprints they leave on the Cosmic Microwave Background (CMB) anisotropies 

Two considerations are in order. Firstly, to explain galactic magnetism it suffices to 
generate "seed" magnetic fields prior to galaxy formation of intensity generally much less 
than 1/iG. In fact, due to magnetohydrodynamic turbulence effects and differential rotation 
of galaxy, extremely week fields can be exponentially amplified. This mechanism, know as 
"galactic dynamo" 0, 8], successfully explains the main characteristics of observed fields 
provided that their intensity and correlation length are B > 10~^'^G and A > lOkpc. If 
dynamo is inefficient, instead, a stronger field, correlated on comoving scales of order IMpc, 
is needed to explain galactic magnetism. In this case the amplification of a primordial seed 
field is just due to magnetic fiux conservation during protogalaxy collapses. Estimates based 
on "spherical infall model" indicate that the time when protogalaxy collapses begins corre- 
sponds to a redshift Zpg not greater than 50 Ij, although galactic disks are assembled at a 
much later epoch, corresponding to redshifts of order z ~ few. Taking the "less conservative" 
value Zpg ~ 50, one can shows that a comoving field B > 10" ^^G explain magnetization of 
galaxies (see, e.g., Ref. IgI). 

Secondly, any viable mechanism of generation must repose on the breaking of confor- 
mal invariance of standard electrodynamics, otherwise the produced fields are vanishingly 
small j^. 

In this paper, we study the possibility to generate seed magnetic fields during infiation 
in a non-invariant conformal theory of electromagnetism described by the Lagrangian 

C = -If,,F^''-^IF,,F'^'', (1) 

where F^^^, = d^A^ — d^A^ is the electromagnetic field strength tensor and F^^ = 
{l/2y/^) e^^pa^Ff"" its dual, with g the determinant of the metric tensor and e^^p^ the 
Levi-Civita tensor. The quantity / is assumed to be a function of a non-trivial background 
pseudoscalar field and its actual form will be specified later. 



In the seminal paper by Turner and Widrow [3|, it was suggested that, during inflation, 
small fluctuating magnetic flelds could be amplifled due to the coupling with the axion fleld 
(in this case, the pseudoscalar J is / = ga-y(t>i where (f) is the axion fleld and Qa-y is the axion- 
photon coupling constant). This possibility was exhaustively studied in the subsequent 
papers by Garretson, Field and Carroll jol, and Field and Carroll [l^. A different scenario 
was analyzed by Prokopec in Ref. [llj where the pseudoscalar was supposed to drive 
(chaotic or extended) inflation. As a matter of fact, in that analysis / was taken to be a 
uniform and slowly varying function of time during inflation and zero subsequently. In all 
cases, however, it was shown that no interesting cosmic flelds can be generated. 

An interaction term of the form IF^j_yF^^ is responsible for the creation of magnetic 
flelds described by an "unbalanced" superposition of left- and right-handed photons. In 
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general, a time variation of / produces a change of F^^F^'^ which can be quantifled by 

j^^d'^xF^yF^" = Hv{t2) — Hv{ti). The quantity Hy is the so-called magnetic helicity 
density (in the volume V) and is proportional to the difference between the number of left- 
and right-handed photons (see Section III). It is deflned by 



V 



d^xA ■ B, 



(2) 



where A and B are the vector potential and the magnetic fleld, respectively. In magneto- 
hydrodynamics [13!], a non- vanishing magnetic helicity indicates a non-trivial conflguration 



of magnetic flux tubes (which are twisted or linked). In particle physics, magnetic helicity 
is known as the Abelian, Euclidean Chern-Simons term. As shown by Jackiw and Pi [l^, a 
helical magnetic fleld results from the projection of a non-Abelian gauge fleld onto a flxed 
direction in isospace. The magnetic helicity coincides then with the winding number carried 
by the non-Abelian vacuum conflguration. 

It is worth noting that, since the magnetic helicity is odd under discrete P and CP 
transformations, if detected on cosmic scales it would indicate a macroscopic P and CP 
violation. This peculiar characteristic of helical magnetic flelds is very attractive and, in 
fact, has induced many authors to devise models of generation of primordial helical flelds in 



the last years 
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mechanisms propos ed, a very intriguing one is that of Cornwall 



also Ref. 



23,^4 



25( 1 . Among the many 



251] ) according to which a net helicity H^ 



V 



and Vachaspati 19l (see 
rib /a PJ| is generated as a by- 



product of baryon-number- violating processes taking place during electroweak baryogenesis. 
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where Ub is the present cosmological baryon number density and a the fine structure constant. 
Interestingly, as it has been pointed out by Giovannini |26|, the reverse it also true: decaying 
non-trivial configuration of (hyper-)magnetic fiux tubes can seed the Baryon Asymmetry of 
the universe. 

The plan of the paper is as follows. In Section II, we consider a toy model which generalize 



the case analyzed by Prokopec in Ref. [ll|] and enables helical seed fields to be produced 
as excitation of the vacuum. In Section III, we compute the amount of magnetic helicity 
associated to the generated fields. In the Conclusions we summarize our results. 

II. CREATION OF HELICAL MAGNETIC FIELDS 
a. General Considerations 

We assume that during infiation the universe is described by a de Sitter metric ds"^ = 
a?{dif' — (ix^), where a{ri) is the expansion parameter, r] = —l/{aH) is the conformal time, 
and H is the Hubble parameter. The conformal time is related to the cosmic time t through 
dt] = dt/a. We normalize the expansion parameter so that at the present time to, a(to) = 1- 

We work in the Coulomb gauge, Aq = Yl^=i ^i^i = 0; '^^ expand the electromagnetic 
field = {Ao, A) as 

A(r7,x)= J2 Aa(r/,x), (3) 

a=l,2 

/d^k 
^k,a ak,aAk,aiv) e'^"" + h.c, (4) 



(27 

where A; = |k| and £k,a (with a = 1,2) are the transverse polarization vectors satisfying the 
completeness relation J2ai^Ka)i{^ta)j — ^ij ~ kikj/k'^. The annihilation and creation oper- 
ators ak,a and aj^^ satisfy the usual commutation relations [ak,o,ak/a'] ~ (27r)^5aQ!'5(k — k'), 
[«k,a, «k',o'] = [c^k a) c^k' a'] ~ 0' with ak,a|0) = 0, whcrc |0) is the vacuum state normalized 
as (010) = 1. 

n 

Introducing the average magnetic and electric fields on a comoving scale A as [5|] 

a^BxAv,^) = Jd'yWx{\^-y\)VxA^{r],y), (5) 
a2EA,«(r/,x)=- jd^yWx{\^-y\)K{v,y), (6) 
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where Wa(|x|) = (27rA^)~^/^e~'^' ^'^'^'^ -* is a gaussian window function and a dot denotes 
differentiation with respect to the conformal time, we can define their vacuum expectation 
values as 

i?L(^) = (0||BU^,x)nO), (7) 

^A,a(^) = (0||EA,«(r/,x)nO). (8) 

Taking into account Eqs. ©-([H]), we obtain 

^Wl{k)VUv), (9) 

-^Wl{k)QUv), (10) 
where Wx{k) = q-^^'^'^Z'^ is the Fourier transform of the window function and 
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-PkAv) = ^,\A,^M\\ (11) 
QkAv) = ^M^M\\ (12) 

are the magnetic and electric power spectra, respectively. 

It is worth noting that in any meaningful model, the function Vk,a/k and Qk,a/k have 
to be integrable in — > (no infrared divergence) in order to have a finite value for the 
magnetic and electric fields -Bf ^, and -Ef q,. 

The equation of motion for A{ri,x.) follows from Lagrangian ([T]): 

A-V^A + /Vx A-V/x A = 0, (13) 



It is worth noting that the effect of the last term in the left-hand-side of Eq. (1131) is just 
to cause a precession of A around V/. Therefore, its presence does not affect the intensity 
of A and, in turn, that of A and that of the related magnetic field. Retaining this term 
would add a major level of complexity to our analysis without, nevertheless, changing the 
final result regarding the average intensity of the inflation-produced magnetic field. For this 
reason, we may write 

A- V^A + /V X A ~ 0, (14) 

where the symbol ~ indicates equality in the sense just discussed. 

We are interested in the study of large-scale electromagnetic fields, that is in modes whose 
physical wavelength is much greater than the Hubble radius, Aphys ^ or equivalently 
\kri\ <^ 1, where Aphys = c^A and A = l//c is the comoving wavelength. 
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If the quantity / is peaked at some small wavelength we expect that, at large scales, the 
term proportional to / in Eq. f|T^ is negligible with respect to the second term. In this case 
we recover the case of free Maxwell theory and no amplification of magnetic modes occurs. 
Therefore, in the following, we consider a simplified special case where I is different from 
zero just at large scales. To see how this assumption modifies Eq. ([Hj), it is convenient to 
work in Fourier space. Inserting Eqs. ([3])-(jl]) in Eq. (fT4|) . we get 



/3=1,2 

where /k is the Fourier transform of / and 

Tk,Q;q,/3 = iq X £q^^ ■ . (16) 

Multiplying Eq. f|T5l) by a\., ^, first on the right and then on the left and subtracting the 
resulting expressions, and integrating in k' and summing on a' afterwards, we obtain 



/3=1,2 ■ 

Since we are assuming that J is different from zero just at large scales ik 0), we may 
write 

/k = (27r)3<5(k)4, (18) 

where (5(k) is a function which is "extremely" peaked at small k. In order to simplify 
calculations, we assume that (5(k) is indeed the Dirac delta function. In this case, only 
modes with small wavenumbers are important when considering the function J^. For this 
reason, we may expand /k-q for small values of k but, due to the presence of the Dirac delta, 
we may expand for small values of q as well: 

/k-q = (27r)='5(k - q)[4 - q • Vk/fc + 0(q')]. (19) 

Inserting the above expression into Eq. (fT7|) . we get to the leading order 

ifc,, + ik^ ± fc4)Afe,„ ^ 0, (20) 

where we used Tk,Q;k,/3 = i^f^a/? and, from now on, ± refer to a = 1, 2, respectively. 

In the next three paragraphs, we will analyze the generation of magnetic fields from the 
vacuum considering a simple phenomenological form for the function Jjt. In Section II e, 
instead, we will give two examples of particle physics models in which that form could 
naturally arise. 
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b. A Simple Toy Model 



Let us suppose that Ik follows a power-law behavior: 

9 



-kif) 



13 



(21) 



where is a dimensionless quantity (which for definiteness we take to be positive) and we 
assume that /5 > 0. [We note that the case analyzed by Prokopec in Ref. corresponds 
to take /5 ^ and gP = const 7^ in Eq. fl^ ]. 

For \kr]\ <^ 1 we can neglect the term proportional to fc^ in Eq. fl2UI) . so that its solution 



IS 



where 



1-P 



(22) 



(23) 



Ci and C2 are constants of integration, and Hi, (x) and Hi, (x) are the Hankel functions of 
first and second kind, respectively. The case {3 = 1 will be analyzed separately (see below). 

Since z oc \kri\^'^~^'^^'^ ^ for |A;?7| ^ 1 we have two different cases: if < /3 < 1 then \z\ <^ 1, 
while if /5 > 1 then \z\ ^ 1. Consequently, in Eq. (l22l) we can use the asymptotic expansion 
of the Hankel functions for small and large arguments, respectively: 



X] 



i2"csc(-Kv) I, I 2 '^[lIfc^cot(7^^^)] u |™| ^ i ^ n 
"T T(l-u) r(l+j/) ' ^ -L, 7^ U, 



< ±2i7r ^Inx, 



2_ g±j(x- 
vra; 



4 2 J 



|x| < l,z/ = 0, 



Ixl > 1. 



(24) 



If < /3 < 1, using the first equation of (j24l) . we find ^^^^(t]) ~ c'l + d2\krj\, where c'^ and 
C2 are constants. Since for modes well inside the horizon {\kri\ 00) we have the plane- 
wave solution 0,(77) = e*'^'' (the normalization corresponds to the standard Bunch-Davies 

j 

vacuum |23]), and for modes well outside the horizon we have the above solution, we can 
fix the values of c'^ and by matching the two solutions and their first derivatives at the 



^ The correct condition in order to neglect the term proportional to k"^ in Eq. ((20|) is: \kri\ <^ {l3gY/'^^^'^\ 
However, neglecting mathematical quibbles such as the case where 5 — > or /3 — > 0, we assume that both 
j3 and g are quantities of order unity, so the above relation reads \kri\ <C 1. 
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horizon crossing, \kri\ = 1. We find c[ = (1 + i)e * and C2 — ^6 ; SO cLS cL final result we 
can write the expression for the electromagnetic field during inflation and at large scales: ^ 

0</3<l: Afc,„(r/) ~ e-^l + ^(1 - M) ]■ (25) 

If /3 > 1, using the third equation of fl24l) . we find Ak^aiv) — |/c?7|'^-'^+^)/^(c'/e*^ + C2e~*^), where 
c'l and C2 are constants. We see that the positive helical modes [corresponding to take the 
plus sign in Eq. ( l23i) ] are ^^^+(77) oc \k'r]\^^^^^/'^ (we neglected inessential oscillating factors), 
while the negative ones are Ak^-{r]) oc l/cr^p+'^^/^e'^'. Therefore, the spectrum of positive 
helicity states is vanishingly small, while that of negative helicity states presents an infrared 
divergence (associated to the exponential factor). Therefore, the case (3 > 1 is meaningless 
and then will be neglected in the following. 

In the case /? = 1, Eq. ( l20l) can be solved exactly (for all \k7]\): 

P = l: AUV) = y|e-(i+2^)/^ M'/'Hi'\M), (26) 

where 

- = ^9 (27) 

and we used the normalization corresponding to the Bunch-Davies vacuum. 

For large scales, we can replace the Hankel function with its small-argument expansion 
in Eq. ([26]). 

c. Electromagnetic Backreaction on Inflation 

Before proceeding further, we want to analyze the problem of backreaction of the inflation- 
produced electromagnetic field on the dynamics of inflation. ^ Since in the following we will 

^ It is worth noting that the matching procedure we have adopted can give us only approximate resuhs. 
Nevertheless, we believe that the simplified analysis performed in the case < /3 < 1 catches the main 
characteristics of the process of creation of magnetic field during inflation. 

^ After inflation, since the conductivity of the cosmic plasma becomes very high and, consequently, the 
electric field is washed out [sj, the electromagnetic energy is simply given by the magnetic one. Moreover, 
the magnetic field evolves adiabatically from the end of inflation until today, and its energy density is 
always subdominant with respect to the energy density of the universe (strictly speaking, this is true for 
cosmic magnetic fields whose intensity is consistent with astrophysical observations, that is i? < 10~^G). 
This means that we can safely neglect the backreaction of the electromagnetic field on the standard 
evolution of the universe from inflation until today. 
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assume that the electromagnetic field does not appreciably perturb the standard evolution 
of the universe, we must verify that the electromagnetic energy density is smaller than 
the energy density associated to infiation which, during de Sitter infiation, is a constant: 
Pinfi = T^- Here, Ti is the so-called reheating temperature, that is the temperature of the 
cosmic plasma at the beginning of the radiation era (here and in the following we assume 
that the reheating phase, during which the energy of the infiaton is converted into ordinary 
matter is "instantaneous" so that, after infiation, the universe enters directly the radiation 
era). 

Starting from the definition of the energy-momentum tensor, 

where S = J d'^x^—g C is the action, and taking into account Lagrangian ([1]), we find 

T,, = \g^.F^pF^^ - F^^F,^, (29) 

that is the standard Mawxell energy-momentum tensor for the electromagnetic field. Here, 
we have assumed, for the sake of simplicity, that the quantity / does not explicitly depend 
on the metric tensor g^^,^ (two particle physics models in which I displays this property are 
discussed in Section II e). The electromagnetic energy density, p = Tq , is then given by 

p = l(E2 + B2), (30) 

where a^B = V x A and a^E = —A. Expanding the electromagnetic field as in Eqs. ([3])-(j4]), 
replacing Bq, and Eq, with their average values Ba,^ and Ea,^, and then taking the vacuum 
expectation value of expression ( l30l) . we obtain the vacuum expectation value of the average 
electromagnetic energy density on a comoving scale A: 



I E i^L + . (31) 



Pa 2 

a=l,2 

In order to have a self-consistent model for the production of a cosmic magnetic field, we 
must verify that, on the scale of interest A, the electromagnetic energy produced during 
infiation is less than the total energy of the universe: 



Pa < Pinfi- 



(32) 
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Taking into account Eqs. (171)- (fT2l) and Eqs. (125|) . (l26!l . and (l3Ti) . we easily find, at large 
scales: 

0<P<1: Pxiv)^^(T) (33) 



^7r2 VA 

and 



p = v. . 4- (i - J |r(.j|^ r (f - (M J"", (34) 



where = y \ + 9- It is important to observe that, in order to avoid an infrared divergence 
in the magnetic power spectrum, we must impose that < 3/2. 

Taking into account that Tf ^ H^mp^ jsl, where mpi ^ 1.22 x lO^^GeV is the Planck 
mass, we have approximatively 



0</5<l: ^ io-i>J ( ) (35) 



-^^10-^°( ) 



and 

P^U^^ 10-'»/(^ ( -I^Y ( M)""; (36) 

where /(x) is an increasing function of x, such that /(1/2) = and /(3/2) ~ 1. Since anal- 
yses of CMB radiation and Big Bang Nucleosynthesis constrain the reheating temperature 
in the range IGeV < Ti < lO^^mpi [3|, and since we are considering the case of large scales, 
\ri\/X <^ 1, we see form Eqs. f l35|) and fl36|) that we can safely neglect the backreaction of the 
electromagnetic field on the dynamics of inflation. 

d. The Actual, Inflation-produced, Magnetic Field 

After inflation, the universe enters the radiation era (assuming instantaneous reheat- 
ing). We restrict our analysis to the case in which, during this era, the interaction term 
in Lagrangian ([1]) is (for some reason) negligible, so that the general expression for the 
electromagnetic field is 

^S(^)=«Me*"' + A,ae-*''- (37) 
Here, ak,a and f3k,a are the so-called Bogoliubov coefficients j^^], determining the spectral 
number distribution of particles produced from the vacuum. By matching expressions (!25|) 
and fl37|) and their first derivatives at the time of the end of infiation, t] = rji, we find the 
spectrum of the electromagnetic field generated from the vacuum at large scales: 

0<(3<1: \Al^M\'=\Pk,a\'^\. (38) 
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In the same way, we find: 



/5 = 1: \AZiVi 



[2--(i-^)r(^)]2 

STrlfcryili+S" ' 



87r|fcr;i| 



1/ > 0, 



iy = 0, 



(39) 



a^+b^cos[2i>ln(|fc7?i|/2)] , c^sin[2t> ln(|fcr;i |/2)] ^ 

47r|fcr)i| 47r|/cr)i| ' * 



Here, > is the imaginary part of u and 



7r 



- (1 + 4z>^)coth(7rz>), 
6^ = (1 - 4z>2)Re[r(z/)2] + 4i>Im[r( 



(40) 



= (1 - 4z>')Im[r(i^)'] - 4z>Re[r(z/)^ 



where Re[x] and Im[x] are the real and imaginary part of x. 

We observe that when g = {u = 1/2) conformal invariance is recovered and then, as it 
should be, we find |A™^(r7i)P = 0. 

From the end of inflation until today, due to the high conductivity of the cosmic plasma. 



■>X,a 



the magnetic field evolves adiabatically [3|, a B 
Inserting Eqs. fl38|) and fl39|) in Eq. ([9]), we obtain respectively 



const, so that B^^"^^ 



< /? < 1 : B 



today 
A, a 



and 



4V27rA2 ' 



(41) 



/3 = 1 : B 



today 



(47r)3/2A2 



/ ln(A/|r;i|) ( \ 



> 0, 



1/2 



(42) 



V2{47r)5/4A2 KM, 

where, for simplicity, we replaced in the third equation of ( !39l) the cosine and sine functions 
with their root-mean-square values 1 / \/2. 
It is useful to observe that 



~ 1.22 X 10^^ — 



A 



mpi 



lOkpc , 



1^1 1 

^ 0.60 X lO-^^A-olp^G, 



(43) 
(44) 
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0.5 1 1.5 2 



g 

FIG. 1: Actual magnetic field in the case /? = 1 as a function of g at the comoving scales A = lOkpc 
(thin lines) and A = IMpc (thick lines), for a reheating temperature Ti = lO^^mpi. Dashed and 
continuous lines refer to the cases a = 1 (positive helicity states) and a = 2 (negative helicity 
states), respectively. The horizontal dotted lines refer to the minimum seed fields required for 
dynamo amplification, B ~ and to directly explain galactic magnetism, B ~ lO^^^G. 



where Aiokpc = A/(10kpc). In Eq. fH5]) . we used the fact that during radiation and matter 
dominated eras the expansion parameter evolves as a oc g^^^^T^^, where T is the tem- 
perature and g*s(T) the number of effectively massless degrees of freedom referring to the 



entropy density of the universe 



28 



Taking into account Eqs. (HTj) and (H3|) - fHi|) . we see that magnetic fields produced in the 
cases < /9 < 1 are well below the minimum seed field required for dynamo amplification, 

B ^ IQ-^'^G. 

In Fig. 1, we show the actual magnetic field in the case = 1 as a function of the 
parameter g at the comoving scales A = lOkpc (thin lines) and A = IMpc (thick lines), for 
Ti = 10~^mpi. Dashed and continuous lines refer to the cases a = 1 (positive helicity states) 
and a = 2 (negative helicity states), respectively. 

As it is clear from the figure, astrophysically interesting fields can be produced as exci- 



4 In this paper, we use the values [281]: Tq ~ 2.35 x lO^i^GeV, ff*s(To) ~ 3.91, and g*s(Ti) = 106.75 
(referring to the number of effectively massless degrees of freedom of Standard Model). It is useful to 
know that IG ~ 6.9 x lO-^OQeV^ and IMpc ~ 1.56 x lO^^GeV^^ 
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tation of the vacuum. They have a definite hehcity ^ and intensity strongly depending on 
the value of g. (The case g > 2, oi equivalently v > 3/2, corresponds to a non-physical 
infrared- divergent magnetic power spectrum). Numerically we find that for g > 0.1 the 
produced field on scales of lOkpc is stronger than the minimum seed field required for a suc- 
cessful galactic dynamo amplification, while if 1.7 < (7 < 2 its intensity on scales of IMpc is 
high enough to directly explain galactic magnetism {B > lO^^^G). In the range of interest, 
0-1 ^ ^ 2, we have approximately 

^toda. ^ ^Q,3.-4r I^^Jl—J Ar-,te G, (45) 



where we used the first equation of fl42l) with a = 2, and Eqs. (I43l) -fl4 



e. Examples of Particle Physics Models 

We want now give two particle physics models in which the function /k takes on the 
form (fTSi) with Ik given by Eq. ( l2Ti) . 

Let us start by considering a model in which I is given by J = 0/M, where M is a mass 
scale and a pseudoscalar (background) field whose Lagrangian density is 

C^ = ^d,<pdy + V{<jy), (46) 

with V{(j)) a potential term that will be specified later. 

In the spirit of mean field theory, we may approximate the function 0(x, t]) by its root- 
mean-square value 4>x{ri) on the scale A: 

0a(^) = (O|| Jd'yW,{\^-yMy,v)\'\0), (47) 

where W^Adxl) is a suitable window function. Introducing the 0-power spectrum, V^{k,ri), 
through the relation 

{OlMvmvm = {27rf5'{k-ci)'^V^{k,v), (48) 
where 0k is the Fourier transform of 0(x), it is straightforward to obtain 

-^m{k)\'V^{k,v), (49) 

^ Since we are considering g > 0, only negative helicity states are astrophysically relevant. Had we taken 
g < 0, the non-vanishing modes would have been the positive ones. 
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where Wx{k) is the Fourier transform of the window function. For the sake of simphcity, let 
us assume that picks out just modes with wavenumber k = 1/X. In other words, 

we take |H/'a(^)P = S{kX — 1), where 6{x) is the Dirac delta function. Then, we may write: 
'PHv) = 'P(l>{^/ ^jV)- Consequently, the Fourier transform of /(x, r^) is simply given by: 

4(,,)^(2.)V(k)IM«!. (50) 

To proceed further, we assume that the </)-power spectrum is peaked at large scales, A — > cxd. 
Taking into account the presence of the Dirac function in Eq. ( !50l) . we can roughly write 
1^(1]) ~ {2n)^5^(k)[P^{k,ri)]^^'^/M, since both k and 1/A go to zero. Therefore, the function 
Ik introduced in Eq. fllSp takes on the form 

To be more concrete, let us now consider the case in which is either a massive pseudoscalar 
field minimally coupled to gravity and a massless pseudoscalar field non-minimally coupled 
to gravity. They are described by Lagrangian fH6|) with 



respectively. Here, is the squared mass of the pseudoscalar field, ^ a real parameter, and 
R the Ricci scalar. It is well-known that in those cases, the 0-power spectrum in de Sitter 
background is given by (see, e.g., Ref. j29|): 

V^{k,v)=(f]\-kvr"'^ (53) 



2n J 



where 



'9 _ m2 

= ^ (54) 
!-12{, 

in the two cases, respectively. Comparing Eq. (!5T|) with Eq. (12T!) . we finally get 

9 = , (55) 

i5 = v^-\. (56) 

It is clear that only negative values of rn? and ^ give positive values of (5. Moreover, the 
astrophysically interesting case discussed in the previous Section, /5 = 1 and 0.1 < < 2, 
corresponds to have rr? = —AH^ and ^ = —1/3, and O.IH < M < l.QH. 
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III. CREATION OF MAGNETIC HELICITY 

Let us now consider the creation of magnetic helicity in our model. The magnetic hehcity 
in a volume = can be conveniently defined as 

= (0| jd^yjdhWx{\^-y\)Wx{\^-z\)A{r],y)-B{r],z)\0), (57) 

where Wa(|x|) is a gaussian window function. Taking into account Eqs. (l3])-(!4j), we obtain 

^Wl{k)n,{v), (58) 

where 

n.iv) = - 1^^,4^)11 (59) 

is the magnetic helicity power spectrum (which is well defined only when Tik/k is integrable 
in A; 0). 

A magnetic field is said to be "maximally helical" at the time f] if either Ak^+iv) or 
Ak~{ff) is zero. 

Looking at Eqs. fl55]) and fISI?]) . we immediately get that in the case < /5 < 1 the 
magnetic helicity is zero. In the case /? = 1, instead, taking into account Eq. (159]) . we obtain 



j^today 



^ (5today)2^ _ (5today)2^^ < g < \, 

(60) 

(Blt^'X - ^ (BX^'^yx, \ <g<l 

where 7± = B(l/2,1 — u) and B{x,y) is the Euler beta function. (The case g > 3/4 
corresponds to an infrared- divergent magnetic helicity power spectrum.) 

It is worth noting that, al least in the case of astrophysical interest, the generated mag- 
netic fields possess a definite helicity (in particular they are left-handed since | <^ 
l^lfc _| ). Hence, they are (almost) maximally helical, and the amount of created magnetic 
helicity is roughly given by 

^today _ _(5today^2^^ (g^^, 

A helical magnetic field leaves peculiar imprints on the CMB radiation. Unfortunately, the 
maximal helicity producible in our mechanism is much smaller than that detectable in near 
future CMB experiments which is, in the most optimistic case, of order of (lO^^G)^Mpc 
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on Mpc scales j30|]. In fact, at those scales we get from Eq. (16T|) and Fig. 1 the value 
~ (10-24G)2Mpc. 



I rrtoday I 



Interesting enough, it has been showed by Kahniashvili and Vachaspati 3l| that, in 
principle, the study of propagation properties of charged cosmic rays through a magnetic 
field could give information about the helicity of the field itself. 



IV. CONCLUSIONS 



Up to today, the hypothesis that all galaxies are pervaded by microgauss magnetic fields 
has not been ruled out yet. Due to their large correlation length, we can entertain the idea 
that they are remnant of an inflationary epoch of the universe. 

In this paper, we have indeed analyzed the production of seed magnetic fields during 
de Sitter inflation. We have considered a photon interaction term in the electromagnetic 
Lagrangian (which breaks conformal invariance of electrodynamics) of the form — ^IF^i,Ff^^, 
where / is a pseudoscalar function. We have considered a simplified special case where / 
is peaked at large scales {k = 1/A 0) during inflation and then is vanishingly small 
afterwards. In particular, we have studied the case where / is parameterized by the power- 
law behavior I{k,r]) = g/{—kri)^, with g a positive dimensionless constant, (3 a real positive 
number, and 1] the conformal time. Also, we have shown that this particular form of / could 
be realized owing to the coupling between the photon and either a (tachyonic) massive 
pseudoscalar field and a massless pseudoscalar field non-minimally coupled to gravity. 

We have found that only when (3 = 1 astrophysically interesting fields can be produced 
as excitation of the vacuum. They are maximally helical and have an intensity strongly 
depending on the value of g (which, in order to avoid an infrared divergence in the magnetic 
power spectrum, is constrained in the range < (7 < 2). In particular, for g ^ 0-1 the 
produced fields on scales of lOkpc are stronger than the minimum seed field required for a 
successful galactic dynamo amplification {B > 10"^'^G), while if 1.7 ^ g < 2 their intensity 
on scales of IMpc is high enough to directly explain galactic magnetism {B > lO^^^G). 



17 



* Electronic address: lleonardo . campanelli@ba. infn. it | 
[1] For reviews on cosmic magnetic fields see: L. M. Widrow, Rev. Mod. Phys. 74, 775 (2002); 

M. Giovannini, Int. J. Mod. Phys. D 13, 391 (2004); D. Grasso and H. R. Rubinstein, Phys. 

Rept. 348, 163 (2001). 
[2] K. Subramanian, PoS MRU, 071 (2007) [arXiv:0802.2"804] [astro-oh]]. 
[3] M. S. Turner and L. M. Widrow, Phys. Rev. D 37, 2743 (1988). 

[4] B. Ratra, Astrophys. J. 391, LI (1992); A. Dolgov, Phys. Rev. D 48, 2499 (1993); F. D. Mazz- 
itelli and F. M. Spedalieri, Phys. Rev. D 52, 6694 (1995); D. Lemoine and M. Lemoine, ibid. 52, 
1955 (1995); M. Gasperini, M. Giovannini, and G. Veneziano, Phys. Rev. Lett. 75, 3796 (1995); 
A. C. Davis and K. Dimopoulos, Phys. Rev. D 55, 7398 (1997); O. Bertolami and D. F. Mota, 
Phys. Lett. B 455, 96 (1999); A. Berera, T. W. Kephart and S. D. Wick, Phys. Rev. D 59, 
043510 (1999); M. Giovannini, ibid. 62, 123505 (2000); 'h ep^ph/0104214| (unpublished); Phys. 
Rev. D 64, 061301 (2001); M. Gasperini, ibid. 63, 047301 (2001); B. A. Bassett, G. Pol- 
lifrone, S. Tsujikawa and F. Viniegra, ibid. 63, 103515 (2001); K. Dimopoulos, T. Prokopec, 
O. Tornkvist, and A. C. Davis, ibid. 65, 063505 (2002); K. Bamba and J. Yokoyama, ibid. 
69, 043507 (2004); ibid. 70, 083508 (2004); A. Ashoorioon and R. B. Mann, ibid. 71, 103509 
(2005); C. G. Tsagas, ibid. 72, 123509 (2005); K. E. Kunze, Phys. Lett. B 623, 1 (2005); 
M. R. Garousi, M. Sami and S. Tsujikawa, ibid. 606, 1 (2005); M. M. Anber and L. Sorbo, 
JCAP 0610, 018 (2006); A. Akhtari-Zavareh, A. Hojati, and B. Mirza, Prog. Theor. Phys. 
117, 803 (2007); K. Bamba and M. Sasaki, JCAP 0702, 030 (2007); K. E. Kunze, Phys. 
Rev. D 77, 023530 (2008); M. Giovannini, Phys. Lett. B 659, 661 (2008); J. Martin and 
J. Yokoyama, JCAP 0801, 025 (2008); K. Bamba and S. D. Odintsov, ibid. 0804, 024 (2008); 
L. Campanelh, P. Cea, G. L. Fogh and L. Tedesco, Phys. Rev. D 77, 123002 (2008). 

[5] T. Prokopec and E. Puchwein, Phys. Rev. D 70, 043004 (2004). 

[6] L. Campanelh, P. Cea, G. L. Fogh, and L. Tedesco, Phys. Rev. D 77, 043001 (2008). 

[7] A. Kosowsky and A. Loeb, Astrophys. J. 469, 1 (1996); J. D. Barrow, P. G. Ferreira, and 
J. Silk, Phys. Rev. Lett. 78, 3610 (1997); R. Durrer, T. Kahniashvih, and A. Yates, Phys. Rev. 
D 58, 123004 (1998); K. Subramanian and J. D. Barrow, Phys. Rev. Lett. 81, 3575 (1998); 
R. Durrer, P. G. Ferreira, and T. Kahniashvih, Phys. Rev. D 61, 043001 (2000); K. Jedamzik, 



18 



V. Katalinic, and A. V. Olinto, Phys. Rev. Lett. 85, 700 (2000); J. Martin and D. J. Schwarz, 
Phys. Rev. D 62, 103520 (2000); J. Martin, A. Riazuelo and D. J. Schwarz, Astrophys. J. 543, 
L99 (2000); T. R. Seshadri and K. Subramanian, Phys. Rev. Lett. 87, 101301 (2001); A. Mack, 
T. Kahniashvih, and A. Kosowsky, Phys. Rev. D 65, 123004 (2002); R. H. Brandenberger and 
J. Martin, Int. J. Mod. Phys. A 17, 3663 (2002); M. Matsumiya, M. Sasaki and J. Yokoyama, 
Phys. Rev. D 65, 083007 (2002); M. Matsumiya, M. Sasaki and J. Yokoyama, JCAP 0302, 003 

(2003) ; L. Campanehi, A. D. Dolgov, M. Giannotti and F. L. Villante, Astrophys. J. 616, 1 

(2004) ; A. Lewis, Phys. Rev. D 70, 043011 (2004); N. Kogo, M. Sasaki and J. Yokoyama, ibid. 
70, 103001 (2004); M. Giovannini, ibid. 71, 021301 (2005); A. Kosowsky, T. Kahniashvih, 
G. Lavrelashvih and B. Ratra, ibid. 71, 043006 (2005); J. Martin and C. Ringeval, JCAP 
0608, 009 (2006); J. D. Barrow, R. Maartens, and C. G. Tsagas, Phys. Rept. 449, 131 (2007); 
M. Giovannini, Phys. Rev. D 70, 123507 (2004); Class. Quant. Grav. 23, 4991 (2006); Phys. 
Rev. D 73, 101302 (2006); ibid. 74, 063002 (2006); PMC Phys. A 1, 5 (2007) [ arXiv:0706.4428] 
[astro-ph]]; Phys. Rev. D 76, 103508 (2007); L. Campanehi, P. Cea and L. Tedesco, Phys. 
Rev. Lett. 97, 131302 (2006) [Erratum-ibid. 97, 209903 (2006)]; Phys. Rev. D 76, 063007 



(2007) ; P. Cea, |astro-ph/070 2293; M. Giovannini and K. E. Kunze, Phys. Rev. D 77, 061301 

(2008) ; ibid. 77, 063003 (2008); ibid. 77, 123001 (2008); D. G. Yamazaki, K. Ichiki, T. Kajino 
and G. J. Mathews, Phys. Rev. D 77, 043005 (2008); F. Finehi, F. Paci and D. Paoletti, ibid. 
78, 023510 (2008); D. Paoletti, F. Finelh and F. Paci. larXiv:0811.0230l [astro-ph]. 

[8] For recent reviews on dynamo mechanisms see, e.g.: A. Brandenburg and K. Subramanian, 



[9 
[10 

[11 
[12 

[13 
[14 

[15: 

[16 

[17: 



Phys. Rept. 417, 1 (2005); A. Shukurov, | astro-ph/0411739| (unpubhshed). 

W. D. Garretson, G. B. Field, and S. M. Carroll, Phys. Rev. D 46, 5346 (1992). 

G. B. Field and S. M. Carroh, Phys. Rev. D 62, 103008 (2000). 



T. Prokopec, [astro-ph/0106247| (unpubhshed). 

L. Campanehi and M. Giannotti, Phys. Rev. D 72, 123001 (2005). 

D. Biskamp, Nonlinear Magnetohydrodynamics (Cambridge University Press, Cambridge, 
England, 1993). 

R. Jackiw and S. Y. Pi, Phys. Rev. D 61, 105015 (2000). 

J. Ahonen, K. Enqvist, and G. Raffelt, Phys. Lett. B 366, 224 (1996). 

J. M. Cornwah, Phys. Rev. D 56, 6146 (1997). 

M. Giovannini and M. E. Shaposhnikov, Phys. Rev. Lett. 80, 22 (1998); Phys. Rev. D 57, 



19 



2186 (1998). 

[18] M. M. Forbes and A. R. Zhitnitsky, Phys. Rev. Lett. 85, 5268 (2000). 

[19] T. Vachaspati, Phys. Rev. Lett. 87, 251302 (2001); Invited talk at International Workshop on 
Particle Physics and the Early Universe (COSMO-01), Rovaniemi, Finland, 30 Aug - 4 Sep 



2001, astro-ph/0111124 



[20] D. S. Lee, W. 1. Lee and K. W. Ng, Phys. Lett. B 542, 1 (2002). 



[21] V. B. Semikoz and D. D. Sokoloff, [ astro-ph/0411496] (unpublished). 
[22] M. Laine, JHEP 0510, 056 (2005). 

[23] A. Diaz-Gil, J. Garcia-Belhdo, M. Garcia Perez and A. Gonzalez-Arroyo, Phys. Rev. Lett. 
100, 241301 (2008). 

[24] C. J. Copi, F. Ferrer, T. Vachaspati and A. Achucarro, Phys. Rev. Lett. 101, 171302 (2008). 
[25] T. Vachaspati, Phil. Trans. Roy. Soc. Lond. A 366, 2915 (2008) [ arXiv: 0802. 15331 [astro-oh]]. 
[26] M. Giovannini, Phys. Rev. D 61, 063004 (2000); ibid. 61, 063502 (2000). 
[27] N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Space (Cambridge University 
Press, New York, 1982). 

[28] E. W. Kolb and M. S. Turner, The Early Universe (Addison- Wesley, Redwood City, California, 
1990). 



[29] A. Riotto, hep-ph/0210162 , Lectures given at ICTP Summer School on Astroparticle Physics 

and Cosmology, Trieste, Italy, 17 Jun - 5 Jul 2002. Published in *Trieste 2002, Astroparticle 

physics and cosmology* 317-413. 
[30] C. Caprini, R. Durrer and T. Kahniashvih, Phys. Rev. D 69, 063006 (2004); T. KahniashviU 

and B. Ratra, Phys. Rev. D 71, 103006 (2005); T. Kahniashvih, Astron. Nachr. 327, 414 

(2006) [ arXiv: astro-ph/05 101 5 1 1 . 
[31] T. Kahniashvih and T. Vachaspati, Phys. Rev. D 73, 063507 (2006). 



